# 定义一个函数，判断一个数是否是素数
def is_prime(n):
  # 如果n小于等于1，返回False
  if n <= 1:
    return False
  # 从2到n的平方根，遍历所有可能的因数
  for i in range(2, int(n**0.5) + 1):
    # 如果n能被i整除，返回False
    if n % i == 0:
      return False
  # 如果没有找到任何因数，返回True
  return True

# 定义一个列表，存储所有的素数
primes = []
# 从1001到2000，遍历所有的数
for n in range(1001, 2001):
  # 如果n是素数，把它加入到列表中
  if is_prime(n):
    primes.append(n)
# 打印列表的长度，即素数的个数
print(len(primes))
# 打印列表中的所有素数
print(primes)
